arXiv: 1501.01135vl [math.CO] 6 Jan 2015 


SOME PROBABILISTIC TREES WITH ALGEBRAIC ROOTS 


OLIVIER BERNARDI AND ALEJANDRO H. MORALES 


Abstract. In this article we consider several probabilistic processes defining random 
grapha. One of these processes appeared recently in connection with a factorization 
problem in the symmetric gronp. For each of the probabilistic processes, we prove 
that the probability for the random graph to be a tree has an extremely simple 
expression, which is independent of most parameters of the problem. This raises 
many open qnestions. 


1. Introduction: an example 

In this paper we consider several probabilistic processes defining a random graph. 
These processes were originally motivated by factorizations problems in the symmetric 
group investigated in [T]. Our main result is a formula, for each of the probabilistic 
processes, of the probability that the random graph is a tree. This probability formula 
turns out to be surprisingly simple, and is in particular independent of most parameters 
of the processes. This is reminiscent of the result of Kenyon and Winkler [3] about 
2-dimensional Branched polymers. We conjecture that more results of this type should 
hold, but we were not able to prove them. 

Before describing our results in full details, let us describe one particular case. Let 
A: be a positive integer. Given a tuple S = {Si, ..., Sk-i) oi k — 1 proper subsets 
of [k] := {1,2,... ,k}, we define G(S) as the digraph having vertex set [k] and arc 
set {ai,... ,afc_i}, where the arc a* has origin i and endpoint the unique integer j in 
[k] \ Si such that {i + l,i + 2,... ,j — 1} C Si, with the integers considered cyclically 
modulo k. For instance, if A: = 7 and S^ = {1,3,6, 7}, then the arc 05 has origin 5 
and endpoint 2. In Figure [T] we have drawn some digraphs G{S) in the case k = 3. 
We now fix a tuple p = {pi,... ,pk) of non-negative integers, and choose uniformly 
at random a tuple S = (iSi,..., Sk-i) of A: — 1 proper subsets of [A:] such that for all 
j e [A:] the integer j is contained in exactly pj of the subsets 5i,..., Sk-i- This gives a 
random digraph G(S). One of the results proved in this paper is that the probability 
that G(S) is a tree (oriented toward the vertex k) is equal to 1 — Pk/{k — 1). This 
result is unexpectedly simple, especially because it does not depend on the parameters 
Pi ?•• • 5 Pk—l • 

Let us investigate in more detail the case A: = 3 of the aforementioned result. The 
situation is represented in Figured} By definition, the tuple S = {Si,S 2 ) is a pair of 
proper subsets of {1,2,3}, and there are three possible trees with vertex set {1,2,3} 
(in general, there are k^~‘^ possible Cayley trees). Let A, B, and C be respectively 
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Figure 1. We represent above the three situations for which the di¬ 
graph G(S) is a tree, with k = 2> and S = {Si, 82 )- 

the events leading to the trees represented in parts (a), (b), and (c) of Figured] For 
instance, 

c = { 2 €Si} n {3 ^ Si} n {3 G S 2 } n {1 ^ S 2 }. 

Now it is not hard to check that the event C has the same probability as the event 
c' = {2 € Si} n {3 ^ ^ 2 } n {3 € Si} n {1 ^ Si}. 

Moreover the events A, B and C' are disjoint and their union is 

AUBUC = ({2 ^ 5i}U{2 € n 3 ^ 5i}U{2 G 5i n 3 G 5i n 1 ^ Si}) n {3 ^ 52 }. 

But since Si is by dehnition a proper subset of {1, 2,3}, the first condition in the above 
clause is always satisfied. It follows that 

F{AUBUC)=F{AUBU C') = P(3 ^ 52 ) = 1 - P 3 / 2 , 

as claimed. Observe that the individual probabilities of the trees represented in Figured) 
do depend on the value of pi and p 2 , but the sum of these probabilities is independent 
of pi and p 2 - 

The rest of the paper is organized as follows. In Section^ we state the main results of 
the paper. In Section [S] we derive a generalization of the matrix-tree theorem tailored 
to our needs: it allows us to express the probability that our random graphs are trees 
as the probability of a “determinant of a matrix of events”. In Section[5]we simplify the 
matrices of events corresponding to our different random processes. In Section [S] we 
compute the determinant of the matrices of events. This computation uses some sign 
reversing involutions. We conclude in Section [6] with some conjectures about another 
random process, and some open questions. 

2. Main results 

In this section we fix some notation and state our main results. We denote by |M| 
the cardinality of a set A. We denote ^ 1 +) i? the disjoint union of two sets A, B. For 
a positive integer k, we denote by [k] the set of integers {1,2, ... ,k}. For i,j G [k], 
we denote by ]i,j] the set of integers (i -|- 1, f -|- 2,... ,j}, where integers are considered 
cyclically modulo k. For instance ]f, i] = 0, ]i, i -|- 1] = {i -|- 1}, and ]i, i — 1] = [A:] \ {i}. 
For an integer r and a tuple p = (pi,... ,Pk) of non-negative integers, we denote by 
Sp^r the set of tuples (S*!,..., Sr) such that for all i G [r]. Si is a subset of [A:] and for 
all j G [A:] the integer j is contained in exactly pj of the subsets Si,... ,Sr- We also 
denote by TZp^r the set of tuples (5i,... , Sr) G such that Si 7 ^ [A:] for all i G [rj. 
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We now define three ways of associating a digraph to a an element in using 
three mappings a, 13, 7 . The mappings a, (3, 7 takes as argument an integer i € [k] and 
a subset S C [k] and are defined as follows: 

• a{i,S) = i if S = [/c], and otherwise a{i,S) is the integer j G [k] such that 
j ^ S but ]i,j — 1] C 5. 

• I3{i, S) = i if S = [k], and otherwise f3{i,S) is the integer j G [k] such that 
j + 1 ^ S' but ]i,j] C S. 

• 7 (z, S) = i if S = [/c], 7(7 S) = i — 1 if f G S, and otherwise 7 (z, S) is the integer 
j G [k] such that j + 1 ^ S but ]i,j] C S. 

The mappings a,j3,'y are represented in Figure [51 




Figure 2. Rules a,j3,^ for creating an arc of the complete graph 
(an informal description of the rule is given between quotation marks 
where i means that i ^ S). 


We now use the mappings a, (3 ,7 to define digraphs. Let S = (Si,..., Sr) G 
and let / be a surjection from [k — 1] to [r]. For Q G {a, j3, 7 }, we define G'^(S, /) to be 
the digraph with vertex set [k\ and arc set A = {oi,..., a^-i} where ai = S/-(j))) 

for all i G [fe — 1]. For instance, the digraph G(S) defined in Section [1] corresponds to 
the case r = k — 1 and G(S) = Ga(S, Id) where Id is the identity mapping from [k — I] 
to [k — 1], Observe that the graph G'(;(S, /) has loops unless S G 7^p,r- 

We are now interested in the probability that the digraph G<^(S, /) is a tree. Observe 
that in this case the tree is oriented toward the vertex k (since every vertex in [k — 1] 
has one outgoing arc). Our main result is the following. 

Theorem 1. Let k and r he positive integers such that r < k, and let p = (pi,... ,pk) be 
a tuple of non-negative integers. Let S = (Si,..., Sr) be a uniformly random element of 
TZp^r (supposing that this set is non-empty), and let f be a uniformly random surjection 
from [k — 1] to [r] independent from S. For f G {a,/3,7}, let P^{p,r) be the probability 
that the random digraph G^(S,/) is a tree. Then 

(a) Pa{p,r) = 1 —pk/r. This is equal to the probability that k ^ Si. 

(h) Pp{p,r) = |7^g,r_l|/|7^p,r| where q = (pi,P 2 - 1,P3 -l,...,Pk- !)• This is 
equal to the probability that Si = {2,3,... ,k}. 
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(c) P^{p ,1 


i=i 


'Pp,r\, where 








[k], = Pi — 1 if i ^ j and q^p = pj. This is equal to the probability that 

\Si\ = k-l. 


Observe that the cardinality of TZp^r appearing in Theorem [T] can be expressed as the 
coefficient of in the polynomial (^0^=1 + 1) “ Y\a=i (^) 

of Theorem [1] was needed to complete the combinatorial proof described in [T] of Jack¬ 
son’s formula [2]. Before embarking on the proof of Theorem[Tl we make a few remarks. 

Remark 1. Theorem [1] can be stated in terms of uniformly random tuples S in iSp^^ 
instead of in 7^p,r- More precisely, for ^ € {a,/3,7}, if S is a uniformly random tuple 
in Sp^r, and / is a uniformly random surjection from \k — 1] to [r] independent from S, 

17? I 

then the probability for the graph G^(S, /) to be a tree is |^^'| r). This is simply 

because the graph f) is never a tree if S € Sp^r \ P-p,r- 

Remark 2. In the case r = k — 1, the result of Theorem [H can be stated without 
referring to a random surjection /. Indeed, for S £ TZp^r and C £ {a, /3, 7 }, let us define 
G^(S) as the graph G(^(S,Id) where Id is the identity mapping from [k — 1] to [k — 1], 
Then the probability T(^(p, k — 1) that the graph G(^(S) is a tree has the same expression 
as in Theorem [TJ For instance, the probability that G„(S) is a tree is 1 — Pk/{k — 1), 
as claimed in Section [T] Indeed, in the particular case r = /c — 1 of Theorem [T| the 
surjection / would be a bijection from [k — 1] to [k — 1] independent from S. But then 
the tuple (S'j(i),..., has the same distribution as S = {Si,..., Sk-i), hence 

the surjection / does not affect probabilities. 

Remark 3. The results in Theorem [1] would hold for any probability distribution on 
the tuples S = (S *!,... ,Sk) of proper subsets of [k], provided that the probability of 
a tuple only depends on the total number of occurrences of each integer i £ [k]. For 
instance, the probability that Ga{S, f) is a tree would still be equal to the probability 
that k ^ Si for such a probability distribution. This result follows from Theorem [1] 
since one can always condition on the total number of occurrences of each integer i £ [k]. 

Remark 4. A slightly weaker version of Theorem [1] can be obtained by not requiring 
the function / to be surjective. More precisely, for ( £ {a,/3,7} and for any positive 
integer r, if S is a uniformly random tuple in TZp^r, and / is a uniformly random func¬ 
tion from [k — 1] to [r] independent from S, then the probability Pd^ip, r) that G,^{S, f) 
has the same expression as in Theorem [TJ For instance, the probability that Gq,(S,/) 
is a tree is l—pkjr. Indeed, this result follows from Theorem [1] by conditioning on the 
cardinality of the image Im(/) of the function /, and by the number of occurrences of 
each integer i £ [k\ in the subsets It is actually this version of Theorem [1] 

(in the case C = 7 ) which was needed in [T]. 
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3. A PROBABILISTIC ANALOGUE OF THE MATRIX-TREE THEOREM 

The matrix-tree theorem is a classical result giving the number of oriented spanning 
trees of a graph as a determinant; see e.g. [1]. In order to prove Theorem [H it is 
tempting to consider the trees on the vertex set [k] as the spanning trees of the complete 
graph Kk, and apply a suitable analogue of the matrix-tree theorem. In this section 
we develop the framework necessary to establish this suitable analogue. 

We first recall the matrix-tree theorem (in its weighted, directed version). Let Kn 
denote the complete digraph having vertex set [n] and arc set A = i ^ j}- We 

call spanning tree of Kn rooted at n a, set of arcs T ^ A not containing any cycle and 
such that every vertex i G [n — 1] is incident to exactly one outgoing arc in T (this is 
equivalent to asking that T is a spanning tree of Kn “oriented toward” the root vertex 
n). We denote by Tn the set of spanning tree of Kn rooted at n. Given some weights 
w{i,j) (taken in a commutative ring) for the arcs (i,j) € A, one defines the weight of 
a tree T G Tn as w{T) = w{i,j). The matrix-tree theorem states that 

{iJKT 

(1) ^ u;(r) = det(L), 

TeTn 

where L = {Lij)i j^[n-i] is the reduced Laplacian matrix, defined by Ljj = —w{i,j) if 
i T j Li^i = u^(l j)- Observe that the matrix-tree theorem gives results 

about the spanning trees of any digraph G with vertex set [n], because one can restrict 
its attention to the spanning trees of G simply by setting w{i,j) = 0 for all arcs (i,j) 
not in G. Observe also that the weights w{n,j) are actually irrelevant in (fT]l. 

We will now derive a generalization of the matrix-tree theorem. Let (0,^, P) be a 
probability space, where H is the sample set, A is the set of events (which is a cr-algebra 
on n), and P is the probability measure. In our applications, 11 will be a finite set and A 
will be the powerset V{Q). Let n be a positive integer and let E = {Eij)i^[n-i],j€[n] be 
a matrix whose entries Eij G A are events. For a tree T G Tn, we define the probability 
of T as 

P(r):=p( fl E.^y 

(*j)eT 

Now we aim at expressing P(T) as some kind of determinant, and this requires 

some notation. Let C[A] be the set of (formal) finite linear combinations of events, with 
coefficients in the held C (i.e., the free C-module with basis A). The elements of C[^] 
are called generalized events and are of the form AiAi A<jAs with Ai,..., A<j G C 

and Ai,As G A. We then dehne the ring R = (C[.4], -|-, n). Here the intersection 
operation “n” is dehned to act distributively on C[Al], that is, for all A,B,G G A, we set 
Ar\{B-GG) = (AnH)-l-(AnC'), and moreover for all A G C we set An(A-H) = A-(AnH). 
We dehne the determinant of a matrix of generalized events M = (Mjj )jjg[„] as 

det(M) = ^ e(7r) n M 2 ,^( 2 ) n • • • n M„,^(„), 

where denotes the set of permutations of [n], and e(7r) is the sign of the permuta¬ 
tion TT. 
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We also extend the probability measure P from A to C[^] by linearity. More con¬ 
cretely, we set P(A-^) = A-P(C') and + B) = P(y4)-|-P(i?). We call ¥-determinant 
of a matrix of generalized events M = and denote by Pdet(M), the prob¬ 

ability of det(M), that is, 

Pdet(M) = P n M 2 ,^( 2 ) n • • • n M„,^(„)) 

Tre&n 

We can now state our generalization of the matrix-tree theorem. For a matrix 
E = {Eij)i^[n-i],j€[n] of generalized events, we define its reduced Laplacian matrix 
L = by setting for all i € [re -1] Lj,j = -Eij ifi^j and Lj,j = 

Proposition 2. Let E = be a matrix of generalized events and let 

L = (Lij)- be its reduced Laplacian matrix. Then 

P(r) =Pdet(L). 

TeTn 

Observe that if the events Eij are all independent, then P(r) = F{Eij), and 

(hj)&T 

Pdet(L) = det {F{Li^j))^ so that Proposition [2] reduces to the usual matrix-tree 

theorem given in ([T|) for the weights w{i,j) = 

Proof. The known combinatorial proofs of the matrix-tree theorem actually extend 
almost verbatim to give Proposition [2l We sketch one such proof, following [5], mainly 
for the reader’s convenience. 

By definition, 

det(L) = Y n 

7rG6„-i ie[n] 

(2) = E n n 

7rG6n-l \i€[n], TT{i)^i j \,iG[n], 7r(i)=i 3 i=-i 

Expanding the right-hand side of ([2]) leads to a sum over a set that we now describe. 
Let Bn be the set of digraphs with vertex set [re] having exactly one outgoing (non¬ 
loop) arc at each vertex i 7 ^ re, and no outgoing arc at vertex re. Let Cn be the set of 
edge-colored digraphs that can be obtained from digraphs in Bn by coloring edges in 
either in blue or red, in such a way that the blue edges form a disjoint union of simple 
directed cycles. We claim that 

det(L) = J] (-i)cycie(c) Pi 

CgCti {i,j)GC 

where cycle(C) is the number of blue cycles of the colored digraph C. Indeed, the blue 
arcs of an element C € Cn encode a permutation vr (the blue arcs are {(z, ' 7 r(z)), for all i 7 ^ 
' 7 r(i)}), the red arcs of C encode a summand in the expansion of n E Ei^j (the 

2€[‘n.],2=7r(2) j^i 

red arcs form a set of the form {{i,j), for all i = TT{i) with j 7 ^ i}), and the factor 

(_l)cycle(C) ig 
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Now, the digraphs in Bn are all the graphs made of a (possibly empty) tree oriented 
toward the vertex n together with a (possibly empty) set of directed cycles on which 
are possibly attached oriented trees. Moreover, if one sums the contribution of all the 
elements C £ Cn corresponding to the same underlying graph B £ Bn one gets 0 if 
there are some directed cycles (because these cycles can be colored either blue or red), 
and otherwise (because all the edges have to be red). This gives, 

det(i) = ^ (_l)«>-d"(c) Pi p Etj. 

CgCti CgTti 

and taking probability on both sides gives Pdet(L) = P(C'). □ 

CdTr, 

Remark. Observe that more generally, for any commutative ring (R,+,0), any 
abelian group ((7,+), and any homomorphism P from {R,+) to ((7,+), there is an 
analogue of the matrix-tree theorem which holds with the same proof: 

^ P( (g) Eij) = P(det(L)). 

TeTn 

Before closing this section we define an equivalence relation on the set C[A] of 
generalized events. Let Ai,..., Am, Bi,..., Bn be events in A, and let Ai,..., Am, 
Pi be complex numbers. We say that the generalized event A = YllLi and 

B = are equivalent, and we denote this A ~ B, if for all w G , the quantities 

= YlT=i and B^ = Y17=i are equal. For instance, for all E,F £ A, 

the generalized events E\J F and E -\- F — E F are equivalent. Also the event 0 and 
the generalized event 0 are equivalent. It is easy to see that ~ is an equivalence relation 
(symmetric, reflexive, transitive) and that \{ E ^ F then P(£') = P(F). Moreover, if 
E F and E' ~ F' then X-E^X-F,E + F ^ E' + F' and E f] F ~ D EX We 
say that two matrices of generalized events E = {Eij)^j^^n] and F = {Fij)ij^[n] are 
equivalent if Eij ~ Fij for all i £ [n\,j G [n]. The preceding properties immediately 
imply the following result. 

Lemma 3. If E and F are equivalent matrices of events, then they have the same 
¥-determinant. 


4. DeTERMINANTAL expressions for the PROBABILITIES P^(p, r). 

We hx r, k, and p as in Theorem [TJ We dehne a probability space (n,M,P) in the 
following way: 

• n is the set of pairs (S, /), where S is in Sp^r and / is a surjection from [k — 1] 
to [r]. 

• M is the power set 'P(fl), 

• P is the uniform distribution on Q. 

We denote by the set of triples (S, /, tt), where (S, /) is in and vr is a permutation 
of [n]. For a matrix of events E = {Eij)ij^[n] we denote by Il{E) the set of triples 
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(S, /, vr) € ^In such that (S, /) is in the intersection nr=i Observe that 

(3) met{E) = ^ Y. 

' ' {sj,7T)en{E) 


We now express the probabilities P^(p,r) defined in Theorem [T] as P-determinants. 
By definition, for € {a,/?, 7 }, P(^{p,r) is the conditional probability, in the space 
(n, A, P), that the random digraph G((S, /) is a tree given that S = (5i,..., Sr) is in 
TZp^r (equivalently, given that none of the subsets Si,... ,Sr is equal to [k ]). Since the 
random digraph G((S, /) is never a tree unless S is in IZp^r (because G^(S, /) has loops 
if S ^ TZp^r) one gets 


^c(p>^) 


P (G^(S, /) is a tree) 
P(S e TZp^r) 


^ J]P(Gc(S,/) = r), 
I I Ter. 


where 7^ is the set of spanning tree of rooted at k. 

We will now use our generalization of the matrix-tree theorem. For in {a,/ 3 , 7 }, 
i € [A; — 1], and j € [fe], we define the event C as the set of pairs (S,/) in O 

such that C(*;5'j(j)) = j. In other words, is the event “the arc a* of the digraph 

G({S,f) is By definition, for any tree T in 7fc, the event G^(S,/) = T is equal 

to r\ii,j)&TE(,i,j. Thus, 


P^(p,r) 



|Pp,r 


E**) n ®<«)' 

TeTk 


Hence by Proposition [2l 

^c(P>^) = ^^IPdet(L^), 
I '^p,r I 


where = (T^^i 7 )jjg[fc_i] is the reduced Laplacian matrix of E(^ = 

We will now define a matrix L'^ equivalent to the reduced Laplacian For t G [r], 
and i,j in [/c] we define the event lE as follows: 


4 = ]i,j]ESt}. 


Observe that lE ■= H becanse 
i A j and jE ■= {(S, /), St = 
generalized events 

T' ^ 


( 4 ) 


T' — 

- 

T' ^ 


]i,i] = 0. We also define the event Jfj by jG = /G if 
[/c]}. For i G [A: — 1] and j G [A;] we define the following 


r/(d _ r/(d 

ij hi-1’ 

t/(*) _ jf{i) 

/■/(d jf{i) r/ld jfE 


Here and in the following, we consider the subscripts of the events lE and jE cyclically 
modulo k] for instance Jqj is nnderstood as J^j. 
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It is easy to check that for in {a,/3,7} and for all i € \k — l\,j G [A:] such that 
i 7 ^ j, one has the equivalence of events ~ instance, 

where 

Ea,i,j = {i + l, j - 1 G and j ^ S’j(i)} = lij-i \ lij, 
hence \ Ii,j ~ hj — h,j-i- Moreover, for all z G [/c — 1], 

Thus, by Lemma [3] the matrices and have the same P-determinant. Our findings 
so far are summarized in the following lemma. 

Lemma 4. For all C G {a,/3,7}, the probability P,^{p,r) defined in Theorem\^is 

PciP^r) = ^^Pdet(L'^), 

where L'^ = (L'^ ^ is the matrix of generalized events defined by (l4ll . 

Next we derive simpler determinantal expressions for the probabilities P^(p,r). 
Proposition 5. For all f G {a,/3,7}, the probability Pi^{p,r) defined in Theorem\^is 

PdP,r) = j^^Pdet(M^), 

where is the matrix of generalized events defined by 

Ma,i,j = Ii!f if and Ma,k,j = 

ifie[k-l], and Mi3^k,j = 

- jfl'lj+i ifi& [k-l], and M^^k,j = 

Example. Let us illustrate Proposition [5] in the case f = a and k = 3. In this case. 



/ 

r/(l) 

^1,1 

r/(l) 

-'1,2 

//(i) 

■'1,3 

\ 

/ n 

2 G 5/(1) 

2, 3 G 5/(1) \ 

II 


r/(2) 

-'2,1 

r/(2) 

-'2,2 

//(2) 

-'2,3 


= 1,3 G 5/(2) 

n 

3 G 5/(2) 


1 

n 

n 

n 

J 

\ n 

n 

n j 


hence by definition of the P-determinant, 

(5) 

Pdet(M„) = P (LI) - P (3 G 5^(2)) + P (2 G 5/(i) O 3 G 5/(2)) - P (2,3 G 5/(i)) 

—P (2 G 5/(1) n 1, 3 G 5/(2)) -I- P (2,3 G 5/(i) PI 1,3 G 5/(2)) • 

Proposition [5] asserts that the probability Pc(p, r) that the graph Ga{S,f) is a tree 
IS I 

is equal to i^pvi Pdet(MQ;). We leave as an exercise to prove that the right-hand side 

17 ? I 

of dS]) is equal to |^^’| x (I — ps/r) as predicted by Theorem [TJ 

The rest of this section is devoted to the proof of Proposition^ We first treat in detail 
the case ( = a. Given Lemma 0] we only need to prove Pdet(L'„) = Pdet(M„), where 
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— ifflij ■ Since P-determinants are alternating in the columns of 

matrices, we can replace the jth column of L'^ by the sum of its j first columns without 
changing the P-determinant. This gives, 

Next, by linearity of the P-determinant in the rows of the matrix, one gets 

Pdet(L'„) = (-l)l^lpdet (M^) , 

DC[fe-l] 

where = (Mf ^ with Mp, = ijP if i G ZD and Mp, = ijP otherwise. We 

now show that only k of the subsets D contribute to the above sum. 


Lemma 6. If D C [k — 1] contains more than one element, then Pdet(M^) = 0. 

Proof. We assume that D contains two distinct integers a and b and want to show 
that Pdet(M^) = 0. We will use the expression (l3|) of P-determinant. By definition, 
a triple (S,/, vr) G rifc-i is in if and only if for all i € D, ]i,k] C and 

for all i G [fe] \ D, ]i, 7 r(i)] C Sf^iy Observe that the above conditions for i = a and 
i = b, namely ]a,k] C Sf(^a) a-iid ]b,k] C Sf(^i,y do not depend on the permutation vr. 
More generally, none of the above conditions is affected by changing the permutation 
TT by TT o (a, b), where (a, b) is the transposition of the integers a and b. Thus a triple 
(S, /, tt) is in Q{M^) if and only if (j){S, f, vr) := (S, /, vr o (a, b)) is in Thus the 

mapping (j) is an involution of Moreover, since the involution (p changes the 

sign of the permutation vr, we get 

Pdet(M^) = -^ Y 

as claimed. □ 

So far we have shown that 

Pdet(L;,) = Pdet(M®) - Y Pdet(M^“>). 

1 ] 

Next, we observe that the set of triples (S,/, vr) in identihes with the set of 

triples (S,/, vr') in n(M„) such that nfk) = k. Indeed, the correspondence is simply 
obtained by replacing the permutation vr of [A: — 1] by the permutation vr' of [k] such 
that TT'{k) = k, and vr'(v) = vr(v) for all z in [/c — 1]. Similarly, for all a G [/c — 1], there is 
a bijection between the set of triples (S, /, vr) in and the set of triples (S, /, vr') 

in such that 7 r'(a) = k. Indeed, the bijection is simply obtained by replacing the 

permutation vr of [fc — 1] by the permutation vr' of [k] such that 7 r'(a) = k, 7r'(k) = 7 r(a) 
and vr'(z) = vr(z) for all z 7 ^ a in [k — 1]. Observe that this bijection changes the sign of 
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the permutation, hence 

Pdet(L;,) = ~ Y 

(S,/,7r)eO(M0) aG[fc-l] (S,/,7r)en(M{“}) 

= Y = Pdet(M„). 

' ' (S,/,7r')Gn(Mc) 

This completes the proof of Proposition [5] in the case C, = a. 

The proof of Proposition [5] in the case ( = P (resp. C = t); is exactly the same as 
the proof given above for ( = a, except that the events ijj are replaced by (resp. 

P — P \ 

5. Computing P-determinants using sign reversing involutions 

In this section we complete the proof of Theorem [T] by computing the P-determinant 
of the matrices Mq for Q € {q;,/3,7}. 

5.1. Computing the P-determinant of the matrix Mq. In this section we compute 
the P-determinant of the matrix M^. 

Proposition 7. The "^-determinant of the matrix M„ defined in Proposition{^ is 

Pdet(M„) = X (l-pk/r). 

Note that Proposition[3together with Proposition [5] prove the case (a) of Theorem[T] 
In order to prove Proposition [TJ we first use the multilinearity of P-determinants with 
respect to rows. For a € [k], we denote by the matrix of generalized 

events represented in Figure [3l that is, = Ma^ij if i > a — 1, = 0ifz<a—1 

and j / i, and if i < a — 1. Note that = Mq. For 

a G {2, 3,... , A: — 1}, we also denote the matrix of events which is 

the same as except the entry in position (a — 1, a — 1) is ^^2;^ (see Figure[3|). 
By linearity of the P-determinant with respect to matrix rows one gets 

Pdet(M(“)) = Pdet(IV(“)) -hPdet(M(“+^)), 

for all a G {2,3,..., A: — 1}. Hence 

k-l 

(6) Pdet(M„) = Pdet(M(2)) = Pdet(M('')) -F ^Pdet(AV(“)). 

a=2 

We will now show that Pdet(M^^^) = x (1 —pk/r)^ and Pdet(Al(“)) = 0 for all 

a G {2,3,... , A: — 1}, thereby proving Proposition [71 

\'JZ I 

Lemma 8. The ^-determinant of the matrix is Pdet(M(^)) = \^~K{l-pk/r). 
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Iff' - Iff' 


rf( 2 ) 


- J; 


/( 2 ) 


0 \ 
0 


V 


0 

rfin) 

ia,l 

r/Ca-l) 

-^a-1,2 

r/(«) 

-^0,2 

n r/(“—2) jf{a—2) 

^a-2,a-2 ’^a-2,a-2 

7-/(a-l) 

^ a — l,a — 2 

^a,a — 2 

0 

E 

j-fM 

-^a,a — l 

jS(a-l) 

jfia) 

J-a,a ’ ' 

0 

^ a — l,k 

jf{o,) 

^a,k 

r/(fc-l) 

^k-1,1 

jfik-l) 

^k-1,2 

Tfik-1) 

^k-l,a-2 

Tf(k-l) 

^k-l,a-l 

jfik-l) 

^k-l,a 

Tf(k-l) 

^k-l,k 


O fl ■■■ fl fi fi ■ ■ ■ fi / 


Figure 3. The matrix (resp. is the matrix represented 

above with the entry E in position (a — 1, a — 1) equal to E = 1^1“ 
(resp. E = 


Proof. Since all the entries in the last row of are equal to 0, one gets 

/k-l 

Pdet(MW) = ^ e(7r)P f| 

Moreover, the only permutations vr contributing to the sum are the identity vr = Id^ 
and the transposition tt = (k — l,k), which gives 



Pdet(M(^)) 


i'k-2 


'k-2 


vi=l 


•* (n (ai? - 4‘”) n j - p I n fiT - -’i 

•" f fi fl? - ■’(?) n (‘PPU - 4'X 


/w _ 7/(d^ 


k—l,k 


V2=l 


By definition, = hi, and is the event = [/c]”. Hence the generalized 

event is equivalent to the event 7 ^ [A:]”. Moreover, 

event “A: G Sf(^k_if\ Hence the generalized event fc-i “ equivalent to 

the event “A: ^ Thus Pdet(M(^)) is equal to the probability of the event “for 

all z G [A: — 1], the subset S'j(j) is a proper subset of [A:] and k ^ . And since / is 

a surjection from [A; — 1] to [r], this event is “the subsets Si,... , Sr are proper subsets 
of [A:] and k ^ S' 7 (fc_i)”. Thus 

Pdet(M(^)) = P(S G TZp^r) X P(A: i S/(fc_i) | S G Up^r)- 

Lastly, the conditional probability P(A: ^ | S G 7^p,r) is equal to 1 —pk/r since 

for any tuple S = (5i ..., Sr) G TZp^r there are exactly pk of the r subsets Si,... ,Sr 
containing the integer k. □ 

It remains to prove that Pdet(A^(“i) = 0 for all a G {2,3,..., A: — 1}. For D C [a — 2], 
we denote by the matrix of generalized events defined by = 
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if i > a - 1 or N^f = if * G D, and N^f = if z G [a - 2] \ D. By 

multilinearity of the P-determinant in the rows of the matrix, one gets 

(7) Pdet(A^(“)) = (-l)l^lpdet(iV“’^). 

riC[a-2] 


It now suffices to prove the following lemma: 

Lemma 9. For all a G {2, 3, ... ,k — 1}, and for all D Cl [a — 2], Pdet(iV“’^) = 0. 

Intuition for the upcoming proof of Lemma[9]. We aim at proving Pdet(A^“’^) = 0 
using ([3|). For this, it suffices to find an involution on the set of triples (S,/, vr) in 
changing the sign of the permutation tt. Now, the hope for finding such 
an involution is based on the observation that the events appearing in the (a — l)th 
and ath rows of the matrix are very similar: for all j in [k] the events N^Li j 

and are respectively “{a}U]a,j] C Sf(^a-i)” and “]a,j] C Sf(^a)\ This implies 

that changing the permutation vr into vr' = vr o (a — l,a) and the surjection / into 
f' = / o (a — l,a) (where (a — l,o) is the transposition of the integers a — 1 and o) 
will only change the requirement “a G Sff^a-i)” iato “a G <S'/(a)”- Now, this can be 
achieved by also transferring the element a from to <5/(0); and this is what 

is done in the proof of Lemma [TT] below. However, this trick does not work for all 
triples (S,/,vr) G because the requirement “a G might still hold for 

(S,/',7r'). Indeed, this happens if there is 6 G [A: — I] such that f{b) = f{a — 1) and 
a g] 6, 7r(6)]. In those cases, one would think about changing vr into vr o (a — 1,6) in¬ 
stead; and this is the intuition behind the proof of Lemmas 1121 and [T3l below. However, 
this does not work in every case, and one has to deal with a few pathological cases 
on the side. We mention lastly that there exist matrices of generalized events having 
a non-zero P-determinant but whose (a —l)th and ath rows coincide with those of 

In order to prove Lemma[9|we need some notation. We hx an integer a G {2, 3,... , k— 
1} and a subset D Cl [a — 2] We dehne &k,a as the set of permutations tt of [/c] such that 
7 r(z) = i for all z G [a — 2]. We also dehne some subsets Kij of [A:] by setting Ki^i = 0 if 
z G [a —2]\D, Ki^i = [k] if z G L)U{a—1}, and Kij =]i,j] if z G {a—1,..., k — l},j G [A:] 
and {i,j) (a — 1, a — 1). By dehnition, a triple (S, /, vr) G Hfc is in if and only 

if the permutation vr is in and for all z G [A: — 1] the subset is contained in 

S'j(j). For TT G / a surjection from [A: — 1] to [r], and t G [r] we dehne the subset 

Rf,iT,t of [A:] by 

-R/,7r,Z = 1^ -ACj 

Observe that a triple (S,/, tt) G Hfc is in if and only if the permutation vr is 

in &k,a and for all t G [r] the subset Rf,TT,t is contained in St- Lastly for j G [A; — 1] we 
denote 

ie[k-l]\{j}, f{i)=f{j) 

SO that Rf, 7 Tj{j) = U Kj .^(^jy 
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We now define a partition of the set by declaring that a triple (S,/, vr) G 

is in 

• W if f{a - 1) = /(a), 

• df if f{a -1) ^ f{a), a ^ and a ^ 

• 3^ if f{a -1) ^ f{a), a ^ and a G 

• Z if f {a - 1) / /(a), a G ilf/,vr,a-i- 

Since = VV’i±)fk’tbl3^l±l^, Equation ([3]) gives 


where 


Pdet(iV“’^) 


W + X + Y + Z 

m 


{ 8 ) W = Y1 *^(^)> Y1 and Z= ^ e(7r). 

(S,/,7r)eW (SJ,7t)&X (SJaIgV (S,/,7r)G2 

We will now show that W = X = y = Z = 0. 


Remark. In the particular case r = k — 1, the surjection / is a bijection and = 0 

for all j. Therefore, in this case X = while W = y = Z = ^. 


Lemma 10. The sum W defined in ([ 8 ]) is equal to 0. 

Proof. We consider the mapping fi defined on W by setting (/)(S,/, tt) = (S,/, vr'), 
where vr' = vr o (a — l,a) (where (a — l,a) is the transposition of the integers a — 1 
and a). Since changes the sign of the permutation tt, it suffices to prove that is an 
involution on W. Clearly this amounts to proving that (S,/, vr') G W. It is clear that 
for all t G [r] \{/(a)}. Moreover, remembering Ka-ij = {a}U]a, j] and 
Ka,j =]a,j] for all j G [k] gives 

^a—l, 7 r(a—l) C ^a, 7 r(a) {®}C ]o, Vr(Q I)]U ]ci, 7r(Q,)] ^a—l, 7 r'(a—l) C -^a, 7 r'(a) i 

hence ftn) = RfTrf(n)- Thus, for all t G H, Rfn't = Rf-nt Y St, that is, 

(S,/,,')£ Sm«/(„ - 1) = /(„) we iJe (S j'.w') e W ae wanted, □ 

Lemma 11. The sum X defined in ([ 8 ]) is equal to 0. 

Proof. We first define a mapping on X. Let (S, /, tt) be in X. Note that a G Sf(^a-i) 
(because a G -fiia-i, 7 r(a-i) ^ 'S'/(a-i))- We now denote by S' = (S'(,..., S’') the tuple 
obtained from S = (Si,..., Sr) by exchanging the presence or absence of the integer a 
between the subsets and Sf(^a)- More precisely. S' = S if a G <S'/(a), and otherwise 

^'f{a) = 'S'/(a)U{a}, = 5/(a_i)\{a}, and S( = Si for all i ^ {/(o-1),/(a)}. We 

now define a mapping on X by setting (j)(S, /, vr) = (S', /', tt'), where S' is defined as 
above, /' = / o (a — 1, a), and vr' = tt o (a — 1, a). We observe the following identities: 

(9) -^a—l,7r'(a—1) ^a,'iT{a) C {o}) R-a,Tr'{a) -^a—l,7r(a—1) \ 

(f-b) R- f'Rf,'K,a, and Hp Rf, 7 T,a—l' 

We will now prove that (j) is an involution on X. Observe that this immediately implies 
that X = 0 because changes the sign of the permutation vr. 
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Claim. The mapping cp is an involution on T. 

Proof of the claim. Clearly, it snffices to prove that if (S',/, vr) is in X then 
{S',f',7r') = (l){S,f,'K) is in X. Let (S,/, tt) G X. Clearly, S' is in Sp^,., f is a 
surjection from [k — 1] to [r], and tt' is in Gk,a- Hence (f{S,f,Tr) is in Next, we 

show that (S',/',7r') is in that is, ^ C S( for all t G [r]. First, for 

t {f{a — l),/(a)}, we have Rf'y,t = Rf,n,t and S( = St, hence Rf/y^t ^ S(. For 
t = f{a — 1) = f'{a), Eqnations ([9| and (fT^ give 

-f^/',7r',a U Ka,Tr'{a) ^f,Tv,a—l U K\ {o} Rf,TT,t \ {n}- 

Moreover, the set S( does not differ from St except maybe for the presence or absence 
of a, hence Rfy^t ^ S(. Lastly, for t = /(o) = /'(a — 1), Eqnations ([9]) and (fT0]l give 

Rf',TT',t Hf ,tt', a—l U Ka—1,1T'{a—l) Hf,TT,a U U {o} U {o}- 

Moreover, S( = StU{a}, hence R^y^t ^ This shows that (S',/',7r') is in 
Enrthermore, (S', /', tt') is in X becanse of (IlOp (and the fact that (S, /, vr) in X). This 

completes the proof of the claim, hence X = 0. □ 

Lemma 12. The sum Y defined in ([5|) is equal to 0. 

Proof. We dehne a partition of y by declaring that a triple (S, /, vr) G T is in 

• Ti if u Kaya) = [k], and u K^yk) = [k], 

• T 2 if Hf,TT,a U Ka,iT{a) ~ [^]) and U -ffa,7r(fc) 7^ [^J) 

Next, we dehne a mapping (/> on T by setting /, tt) = (S, /, vr'), with 

• tt' = TT o {a, k) if (S, /, tt) G Ti, 

• tt' = TT o (6, k) if (S, /, tt) G T 2 , 

• tt' = TT o (a, b) if (S, /, tt) G Ta, 

where b is the greatest integer in [A; — 1] \ {a} snch that f{b) = f{a) and a G 
(snch an integer b necessarily exists since a G Hf^T^^a)- Since the mapping clearly 
changes the sign of the permntation tt, showing that cj) is an involntion on y wonld 
imply Y = 0. In the rest of this proof we show that (j) is an involntion on y and in fact 
an involntion on each of the snbsets Ti, T 2 , 3 ^ 3 - 

Claim 1. The mapping described above is an involntion on Ti- 

Proof of Claim 1. It clearly snffices to prove that if (S, /, tt) is in Ti then (S, f, tt') = 
(p{S, f, tt) is in Ti- Eirst observe that by dehnition of Ti, 

Rf,nJ(a) = m = Rfyj(a)- 

Hence Rf^TT,t = Rf,ir',t for all t G [r], Thns Rfy,t = Rf,iT,t ^ St, that is, (S,/, vr') 
is in H(A^“’'^). Next, we show that (S,/, vr') G T- Since (S,/, tt) is in y we get 
a Hf,TT,a = Hfy,a-i- Moreover, a G Hfy^a becanse a G [k] = Rfyjy but 
a i Kay{a)- Thns (S,/,7r') G y. Lastly, (S,/,7r’) G Ti clearly implies (S,/,7r') G Ti- 
This completes the proof of Claim 1. 
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Claim 2. The mapping (p is an involution on 3 ^ 2 - 

Proof of Claim 2. Let (S, /, vr) be in T 2 - We need to prove that (S, /, vr') = 4>{S, f, vr) 
is in 3^2 and that b is also the greatest integer in [L — 1 ] \ {a} such that f{b) = f{a) 
and a G First observe that b > a — 2, otherwise we would have Tr{b) = b and 

^b,TT{b) equals 0 or [k] but either is impossible (indeed 7 ^ 0 because a G 

and / [k] because (S, /, vr) is in T 2 ); hence in particular vr' = vro ( 6 , k) is in &a,k- 

We now show that (S,/, vr') is in that is, Rfy^t ^ St for all t G [r]. First, 

for t 7 ^ /(a), one has Rfy,t = Rf, 7 r,t ^ St- Second, for t = f{a), we have 

[A;] Hf,n,a U Ka,TT{a) Rf,TT,f{a) — 

Hence, Sj^a) = [k] and Rfyj(a)^ C Sfy Thus (S,/, 7 r') is in 

We now prove that (S, /, vr') is in T 2 - Let c be the unique integer in [k] \ such 

that ]c, a] C (tLe integer c necessarily exists since a G Hf,TT,a 7 ^ [fc])- We will 

show 

(11) a < 7r{k) < c < b < k, and a < 7 r( 6 ) < c. 

Since b ^ [a — 2], and b ^ {a — l,a}, we get a < b < k. Moreover, by defini¬ 

tion of b and c, we get ]b,a] C]c,a], hence a < c < b. Furthermore, we know 
]c, a]U ]a, 7 r(/c)] C Lfj ]a, 7 r(/c)] 7 ^ [k], hence a < ’K{k) < c. Lastly, we know 
a g] 6 , 7 r( 6 )] and ]c, 6 ]U] 6 , 7 r( 6 )] C ^ [k], hence a < 7 r( 6 ) < c. Thus (flTTl holds. 

It implies a g] 6 , 7 r(/c)] C Hfy^a- Thus (S,/,vr') is in T. Moreover c ^]b,7r{k)], hence 
c ^ Hfy^a- Thus Lfp 7 r',a 7 ^ [k]- We also know ]«, 7 r(a)] = [A:], hence c G ]a, 7 r(a)], 

which implies ]a,c] C]a, 7 r(a)]. Moreover, ]c,b] C]c,a] C and 7 r( 6 ) < c, hence 

]c,b]C [J Kiyty Renee, 

ie[A:-l]\{a, 6 }, f[i)=f{a) 

L^/,7r',aU]a,7r(a)] D]c,b]U]b,a]U]a,c] = [k]. 

Thus (S, /, it') is in 3 ^ 2 - Lastly, we know a G ]b, vr(A:)] = Kijy(iy, hence b is the greatest 
integer in [A: — 1] \ {a} such that f{b) = /(a) and a G This shows that p is an 

involution on T 2 and completes the proof of Claim 2. 

Claim 3. The mapping p is an involution on T 3 . 

Proof of Claim 3. Let (S, /, vr) be in T 3 . We need to prove that (S, /, vr') = p{S, /, vr) 
is in T’s and that b is the greatest integer in [A; — 1] \ {a} such that f{b) = f{a) and 
a G First observe that 6 ^ [a — 2] (indeed 7 ^ 0 because a G and 

Rb,iT{b) 7 ^ [k] because (S, /, vr) is in T 3 ), hence vr' = vr o (a, b) is in 6a,k- We now show 

(12) a < 7 r(a) < b < k, and a < Tr{b) < b. 

Since b ^ [a — 2], and 6 7 ^ {a — 1, a}, we get a < b < k. Moreover, we know a g] 6 , 71 ( 6 )], 
hence a < 7r{b) < b. Moreover, we know ]a, 7 r(a)]U] 6 , 7 r( 6 )] C LApjr,aU]a, 7 r(a)] 7 ^ [A:], 
hence a < 71 ( 0 ) < b. Thus (1121) holds. It implies 

(13) ]a, 7 r'(a)]U ]6,71^(6)] = ] 6 , max( 7 r(a), 7 r( 5 ))] = ]a, 7 r(a)]U ] 6 , 7 r( 6 )] 

hence Rfyjy = Rf,-K,f{a)- Thus, for all t G [r], Rfy^t = Rf, 7 T,t C St, that is, (S,/, 7 r') 
is in H(AV“’-^). 
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Equation (fT^ also implies a €] 6 , 7r(a)] C Hfy^a- Thus (S,/, vr') is in y. Moreover, 
([ni) gives Hfy^a^]a,TT'{a)] = ]a, 7r(a)] ^ [k]. Thus (S,/, tt') is in Ts- Lastly, 

a G ]b, 7r(a)] = hence b is the greatest integer in [/c— l]\{a} such that f{b) = f{a) 

and o € Kby{b)- This shows that (j) is an involution on Ts and completes the proof of 
Claim 3. 

Claims 1,2,3 imply Lemma [T 2 j □ 

Lemma 13. The sum Z defined in ([ 8 ]) is equal to 0. 

Proof. The proof of Lemma [13] is very similar to the proof of Lemma [12] (actually, it is 
identical except “a” is replaced by “a — 1” in certain places). We first dehne a partition 
of Z by declaring that a triple (S, /, vr) G Z is in 

• -^1 T U [A;], and U Ka—i,'K{k) [^ji 

• Z 2 if 1 U = [A;], and U Ka—i,n{k) 7 ^ \fi\i 

• if 1 U ^ [A;]. 

Next, we define a mapping on Z hy setting ^(S, /, vr) = (S, /, vr'), with 

• tt' = TT o (a — 1, A:) if (S, /, vr) G Zi, 

• tt' = TT o ( 6 , k) if (S, /, tt) G Z 2 , 

• tt' = TT o (a - 1, 6 ) if (S, /, tt) G ^ 3 , 

where b is the greatest integer in [A: —l]\{a—1} such that f{b) = /(a—1) and a G 
(such an integer b necessarily exists since a G LAj^ 7 r,a-i)- Since the mapping clearly 
changes the sign of the permutation tt, showing that is an involution on Z would 
imply Z = 0. Actually, the proof that is an involution on each of the subsets Zi, 
2 ^ 2 ) and - 2 I 3 is almost identical to the one of Lemma [12] (for instance, the identities (llip 
and (fT^ still hold with c the integer in [A:] \ LA/_ 7 r,a-i such that ]c, a] C fAj ,^, 0 - 1 ) and is 
left to the reader. □ 

Lemmas [H] [121 [E] imply Lemma [9] Thus, by ([7]), Pdet(A^*^“)) = 0 for all a G 
{2, 3,... , A: — 1}. This together with Equation ([6]) and Lemma [8] complete the proof of 
Proposition [7] hence of the case (a) of Theorem [1] 

5.2. Computing the P-determinant of the matrix Mg. In this section we compute 
the P-determinant of the matrix M^. 

Proposition 14. The ¥-determinant of the matrix defined in Proposition is 

Pdet(M^) = 

where q= {pi,P 2 - 1 ,P 3 - 1,... - !)■ 

Note that Proposition [TT] together with Proposition [5] prove the case (6) of Theo¬ 
rem [1] We now sketch the proof of Proposition [TT] which is very similar to the proof of 
Proposition [7] We first use some column and row operations on the matrix M^. Let 
where M^if f / A: and M^ ;,^ = Q. Since the matrix 
M'p is obtained from by reordering its columns cyclically we get 

Pdet(M^) = (-l)^-^Pdet(M^). 



18 


OLIVIER BERNARDI AND ALEJANDRO H. MORALES 


( ■ 

jfw 

'^l,a 

t/(1) 

"^l.a + l 


7-/(a-l) 

"^<1 — 1,1 
jfG) 

'Ja,l 

0 

E 

0 

jfG-U 
■^a-l.a + l 

7 /(a- + l) rf{a+l) p. 

*^a,+ l,a+l •^a,+ l,a+l ^ 

jf(a) 

•^a.k 

0 


V 


0 

n 




Q 


T/(fe-l) _ Tf(k-l) 

o 


0 

o 




Figure 4. The matrix (resp. is the matrix represented 

above with the entry E in position (a, o) equal to E = (resp. 


E = 


For a G [A: — 1], we denote by [fc] the matrix of generalized events 

represented in Figure 01 that is, ■ j = if i < a, = 0 if i G 

{a + 1,..., fe - 1} and j / i, and if i G {a + 1,..., A; - 1}, 

and = Q. Note that = M^. For a G {2, 3,... , A: — 1}, we also denote 

the matrix of events which is the same as except the entry in 
position (a,a) is (see Figure 0|). 

By linearity of the P-determinant with respect to matrix rows one gets for all a G 
{2,3,...,A:-1}, 

Pdet(M(“)) = Pdet(AV(“)) +Pdet(M(“-^)). 

Thus, using this relation iteratively starting with Pdet(M^) = Pdet(M(^“^)) we get 

fc-i 

(14) Pdet(M^) = Pdet(M(^-i)) = Pdet(M(i)) + ^ Pdet(iV(“)). 

a=2 

We will now prove that Pdet(M^^^) = (—p, and Pdet(A^*^“)) = 0 for all 

|<^p,r I 

a G {2,3,... , A; — 1}. 


Lemma 15. The F-determinant of the matrix appearing in m is Pdet(M(^)) = 

( pfc-i 

Proof. This proof is similar to the one of Lemma El One easily gets 

^k-l 


Pdet(MW) = F 

\i=2 / 

= (-1)*->P (n {l!f> - j(P) n (j{W - j/W)) . 


\i=2 
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Moreover, the generalized event inside the argument of P(-) is clearly equivalent to the 
event = {2, 3,... , A:} and for all z G {2,..., A: —2}, S'j(j) is a proper subset of [A;]”. 

This event is the same as “<5/(1) = {2,3,..., A:} and S € 7^p,r” which has probability 

P(S G n^,r) X nSfil) = {2,3,..., A:} I S G 7^p,,) = x 

|^p,r| |'^p,r| |^p,r| 

□ 

It remains to prove that Pdet(A^^“)) = 0 for all a G {2,3,... , A: — 1}. The proof is 
very similar to the one presented in Section 15.11 (indeed we point out that the rows a — 1 
and a of are identical to the ones in Section [5. ip . For DC {a + 1,..., A: — 1}, we 
denote by the matrix of generalized events defined by if 

i < a OT j i, and if i € D, and if z G {a + 1,..., A: — 1} \ D. 

By linearity of the P-determinant, one gets 

Pdet(lV(“)) = (-l)l^lpdet(iV“’^). 

DC{a+l,...,/c—1} 

We now fix a G {2,3,..., A: — 1} and D C {a + l,...,A: — 1} and proceed to prove that 
Pdet(A^“’^) = 0. We define as the set of permutations vr of [A:] such that 7r(z) = z for 

all i > a. We also define some subsets Kij of [A:] by setting iC/j = [A:] if z G D U {a — 1}, 

Ki^i = 0 if z G {o, o + 1,..., A: — 1} \ D, and Kij = ]i,j] for i G {1, 2,, a}, and 

j ^ [^] \ {*}■ By definition, a triple (S,/, tt) G 12^ is in if and only if the 

permutation tt is in ©^, and for all z G [A: — 1] the subset is contained in 5'/(j). 

Lastly for / G [A: — 1] we denote 

We now define a partition of the set f2(A^“’^), by declaring that a triple (S,/, vr) G 
is in 

• W if f{a - 1) = /(a), 

• T if /(a - 1) 7 ^ /(a), a ^ and a ^ 

• T if /(a - 1) / f{a), a ^ and a G LA/, 7 r,a, 

• Z if /(a - 1) / /(a), a G 

Since = Vyi±)T’l±)T'WZ, we get Pdet(A^“’^) = > where 

W= ^ e(7r),Y= and Z= ^ e(7r). 

{S,/,7r)GW {S,/,7r)GA (S./,7r)e3^ (S,f,7r)eZ 

We then show that W = X = Y = Z = 0. 

The proof that IT = 0 is identical to the proof of Lemma [lOl The proof that X = 0 
is identical to the proof of Lemma fTTl It is done by defining a sign reversing involution 
(j) on X. Explicitly, the mapping (j) is defined on X by setting ^(S,/, vr) = (S',/',7r'), 
where /' = / o (a — 1, a), and vr' = tt o (a — 1, a), and is obtained from S simply by 
exchanging the presence or absence of the integer a between the subsets Sf(^a-i) and 
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Sf(^a)- The proof that (p is an involution on X is identical to the one given in the proof 
of Lemma [TT] 

The proof that T = 0 is identical to the proof of Lemma [T2j It is done by defining a 
sign reversing involution (p on y. Explicitly, one considers the partition y = 
dehned exactly as in the proof of Lemma [T2j One then defines a mapping p on y by 
setting p{S, /, vr) = (S, /, tt'), where tt' = tt o (a, k) (resp. tt' = tt o (6, k), tt' = vr o (a, b)) 
if (S, /, tt) is in Ti (resp. ^ 2 ; Ts), where b is the greatest integer in [/c — 1] \{a} such that 
f{b) = f{a) and o G iLfe ,r(6)- The fact that p is an involution on y (actually, on each 
of the subsets Ti, 3^2, and Ts) is identical to the one given in the proof of Lemma fT^ 
The proof that Z = 0 is again identical to the proof of Lemma [131 This completes 
the proof that Pdet(A^“’^) = 0 for all a G {2,3,..., A: — 1} and all D C {o+l,..., A: — 1}. 
This together with LemmafTHl complete the proof of Proposition fTdl hence of the case (b) 
of Theorem [H 

5.3. Computing the P-determinant of the matrix M^. In this section we compute 
the P-determinant of the matrix M^. 

Lemma 16. The ¥-determinant of the matrix defined in Proposition\^is given by 

k 

Pdet(M.^) = ^Pdet(Q(“)), 

a=l 

where matrix of generalized events defined by if 

i G [a - 1], ifie{a + l,...,k}, and = 0. 

Note that = Mg, where is the matrix defined in Proposition [5l 
Corollary 17. The ¥-determinant of the matrix defined in Proposition\^is 

Pdet(M.,) = ^^f^^, 

where ..., and for all i G [A:], = pi — 1 if i j and = pa- 

Proof of Corollary By Proposition [TH 

Pdet(Q(*^)) = Pdet(M^) = 

Now for a G [A; — 1], permuting cyclically k — a times both the rows and the columns 
of the matrix <5^“^ gives 

Pdet(Q(“)) = Pdet(i?(“)), 

where is the matrix defined by for i G [A: — a], 

^i,j ~ -^i+i^^+a+i for i G {A: — a + 1,... , A: — 1}, and = 11. Since the surjection / 
is uniformly random we get. 


Pdet (!?(“)) = Pdet(T(“)) 
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where is the matrix defined by for i € [k — 1] and 

Tjfj = Q. Now, by symmetry (obtained by replacing the integer i by i + a in the subsets 

Si,..., Sr), the P-determinant Pdet(T(“)) is obtained from Pdet(Ma) by replacing pi 
by pi^a for all i & [k], that is, 

Pdet(r(“)) = 

This together with Lemma [TBl completes the proof of Corollary 1171 □ 


Note that Corollary 1 171 together with Proposition [5] proye the case (c) of Theorem [T] 
It now only remains to proye Lemma fTHl For DC [fc — 1] and j G [k], we denote = 

the matrix of eyents defined by MP- = if i G H, MP = if 


mp 

hj 


i,je[k] 


i ^ [k — 1]\D, and = 11. By the multilinearity of P-determinants with respect to 
the rows, 

Pdet(M^) = (-l)I^IPdet(M^) . 

DC[fc-l] 

We now show that only k of the subsets D contribute to the aboye sum. 

Claim. If H C [/c — 1] contains an integer a > 1 but not a — 1, then Pdet (M^) = 0. 


Proof of the claim. We suppose that the subset D contains the integer a > 1 but 
not a — 1. We remark that the (o — l)th and ath rows of the matrix are almost 
identical: indeed for all j G [fe], ■ = jlYj+i while M^- = We now 

show that this implies Pdet(M^) = 0. We define a mapping 0 on by setting 

(^(S, /, vr) = (S, / o (o — 1, a), vr o (a — 1, a)). We want to proye that (p is an inyolution 
on Q{M^). Let (S, /, tt) G Hfc. It is clear that / o (a, a — 1) is a surjection from [k — 1] 
to [r] so that (p{S,f,TT) is in Moreoyer, it is easy to see from the aboye remark 
about the (o — l)th and ath rows of that (t){S, f, tt) is in Q{M^) (indeed the triple 
4>{S,f,7r) satisfies the condition because (S,/, vr) satisfies and 

(p{S,f,7r) satisfies because (S,/, vr) satisfies Thus the mapping 

(j) is an inyolution of ll(M^). Moreoyer, since the inyolution p changes the sign of the 
permutation vr we get 

Pdet(M^) = -^ Y 

(S,/,7r)GlI(MO) 


as claimed. 


So far we haye proyed that 


k-l 

Pdet(M.^) = Pdet(M®) + ^(-1)“ Pdet(M[“l). 

a=l 
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IVfnrpnvpr Liv rlpfinitinn Af® = anrl it, is pasv tn spp hv rpnrrlprino' ttip rnws nf 



the proof of Lemma [TUI hence of the case (c) of Theorem [TJ 

6 . Open questions and conjectures 

This paper leaves open a few questions. First of all, the proof we obtained of Theo¬ 
rem [1] is not as elegant as we had hoped for, hence the following question: 

Question 18. Is there a more direct proof of TheoremUl? In particular, can one prove 
Theorem{^ without using the matrix-tree theorem? 

Second, our use of the matrix-tree theorem (for which a “forest version” exists) 
suggests the following question: 

Question 19. For C € {a, (3, 5}, is there a simple expression for the probability that the 
random graph G({S, f) is a pseudo-forest, that is, a digraph in which the only cycles 
are loops (so that the connected components are trees oriented toward a vertex which 
has a loop)? 

More importantly. Theorem [1] gives simple formulas for the probability that certain 
random graphs are trees. We stumbled over one of these formulas while studying a 
factorization problem in the symmetric group [T], but it is unclear how general this 
phenomenon is. 

Question 20. Is there a more general theory encompassing Theorem[l\? 

We conclude with a few conjectures about a random digraph Gs defined analogously 
to Ga, Gjs, and G^. Let 6 be the mapping taking as argument an integer i G [k] and a 
subset S C [A:] and defined by 6{i,S) = i i G S, and otherwise 6{i,S) is the integer 
j G [/c] \ S' such that ]i,j — 1] C S. The mapping 6 is represented in Figure El 



kt € 5 or 5 = [A:] \ {i} 


i + l € S, i,i+2 ^ S 


(5) Rule • 1 


Figure 5. Rule 5 for creating an arc of the complete graph Kk- 


Given a tuple S = (Si,..., Sr) in Sp^r and a surjection / from [k — 1] to [r], we 
define Gs{S,f) to be the digraph with vertex set [k] and arc set A = {ai,... , 0 ^- 1 } 
where a* = (i, (5(i, Sj-(j))) for all i G [A: — 1]. The following conjecture suggests there is 
a simple expression for the probability that ^^(S, /) is a tree. 

Conjecture 21. Let k and r be positive integers such that r < k, and let p = 
{pi,... ,Pk) be a tuple of non-negative integers. Let S = (Si,..., Sr) be a uniformly 
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random element oflZp^r (supposing that this set is non-empty), and let f be a uniformly 
random surjection from [k — 1] to [r] independent from S. Then the probability Ps{p, r) 
that the random digraph Gs{S,f) is a tree is 

Ps[P,r) = —^ 

I '^p,r I 

This is equal to the probability that S'! = 0. 

Although we were unable to prove ConjectureEH we obtained the following analogue 
of Proposition [5l 

Ps{p,r) = ^E£^Pdet(M,), 

I I 

where Mg = is the matrix of the generalized event defined by = 

j for i G [k — 1], and = Q. Thus Conjecture ED is equivalent to the 

following identity: 

(15) Pdet(M5) = 

|Op,r| 

One approach to try to prove (1151) is to use the multilinearity of P-determinants. For 
D C [k — 1] and j G [k], we denote = (mP^ , where MP = ii i G D, 

MPj = if i G [A: — 1] \ H, and = Tl. By multilinearity of P-determinants, 

(16) Pdet(M5) = (-l)l^lpdet (M^) . 

DClfc-l] 

Calculations suggest the following formula for the P-determinant of M^-. 

Conjecture 22. For all DP [A: — 1], 

Pdet(M^) = x¥{k ^ Si,D C Si \ S G 7^p,^). 

Remark. Conjecture [22] implies Conjecture 1211 since substituting the conjectured for¬ 
mula for Pdet(M^) in (|16p and doing inclusion-exclusion gives 

nMs) = X F{Si = 0 I S G 7^p,,). 

Remark. When H = 0 then = Mq, as defined in Proposition [5] Hence, in this 
case Conjecture [22] holds by Proposition [7] Similarly, when D = [A: — 1] the matrix 
is the matrix Mjg as defined in Proposition [5] up to shifting cyclically the k columns 
one position to the right and shifting the first A: — 1 rows one position down. Hence, in 
this case Conjecture [22] holds by Proposition [TT] 
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